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Abstract 

We present the results of a search of log-periodic corrections to scaling in the moments 
of the energy dissipation rate in experiments at high Reynolds number (w 2500) of three- 
dimensional fully developed turbulence. A simple dynamical representation of the Richardson- 
Kolmogorov cartoon of a cascade shows that standard averaging techniques erase by their very 
construction the possible existence of log-periodic corrections to scaling associated with a dis- 
crete hierarchy. To remedy this drawback, we introduce a novel "canonical" averaging that 
we test extensively on synthetic examples constructed to mimick the interplay between a weak 
log-periodic component and rather strong multiplicative and phase noises. Our extensive tests 
confirm the remarkable observation of statistically significant log-periodic corrections to scal- 
ing, with a prefered scaling ratio for length scales compatible with the value 7 = 2. A strong 
confirmation of this result is provided by the identification of up to 5 harmonics of the funda- 
mental log-periodic undulations, associated with up to 5 levels of the underlying hierarchical 
dynamical structure. A natural interpretation of our results is that the Richardson-Kolmogorov 
mental picture of a cascade becomes a realistic description if one allows for intermittent births 
and deaths of discrete cascades at varying scales. 

1 Introduction 

Since Kolmogorov's presentation of his theory of three-dimensional fully developed turbulence in 



1941 (see [13] for a modern account), there has been a vigorous continuing investigation aimed 
at verifying its accuracy and testing for possible deviations. While its main predictions (/c -5 / 3 
spectrum of velocity fluctuations, universality of the exponents of structure functions in the inertial 
range) are verified to a good approximation, there are clear deviations. The most notable is the 
nonlinear dependence of the exponents Q q of the structure functions as a function of their order q, 
which is a hallmark of multifractality [[35|]. There is also a strong interest in testing Kolmogorov's 
universality assumptions that all the small-scale statistical properties are uniquely and universally 



determined by the scale £, the mean energy dissipation rate e and the viscosity (see [38] for recent 
evidence of the contrary) and that there is a well-defined limit at infinite Reynolds number (see 
ITCL ff8|] for alternative scenarios). 



Kolmogorov's K41 hypotheses "were bas ed physically on Richardson's idea [ |39| ] of the exis- 



tence in the turbulence flow of vortices of all possible scales..." [27] and turbulence phenomenology 



is often cast in terms of a cartoon depicting a cascade of the energy introduced into the largest eddies 
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of size lo and "cascading" down a hierarchy of eddies of size £ n = io/^y 71 where 7 > 1 (often taken 
equal to 2 but without particular significance: [13] bottom of page 103) at the same rate e and being 
eventually removed by dissipation at the smallest scales. To our knowledge, Novikov was the first 
to point out in 1966 that, if the discrete hierarchy of the cascade is taken seriously with a specific 7, 
the structure functions in turbulence should contain log -periodic oscillations [[33], fPj ]. His argument 
was that, if an unstable eddy in a turbulent flow typically breaks up into smaller eddies of about half 
or a third of its size (7 = 2 or 3), but not into eddies ten or twenty times smaller, then one can 
suspect the existence of a preferred scale factor 7, hence the log-periodic oscillations. Smith et al. 
[43] confirmed on explicit geometrical examples involving generalizations of the triadic Cantor set 
that log -periodic corrections to scaling arise in fractal models with lacunarity, i.e., having a pref- 
ered scaling ratio 7 between scales. Shell models construct explicitly a discrete scale invariant set 
of equations on a hierarchy of shells k n = 2 n ko in momentum space whose solutions are marred 
by unwanted log-periodicities (see for instance [|1|] for a refinement addressing the problem). 

Oscillations in log-log plots of structure functions and other quantities are often observed but 
do not seem to be stable and depend on the nature of the global geometry of the flow and recir- 



culation [ |13 , g] as well as the analyzing procedure. In his book ([]13[], pages 130-131), Frisch 
discusses the possible existence of such log-periodic corrections and points out that the "extended 
self-similarity" technique ^ which has significantly improved the accuracy of exponents may be 
useful in particular because it overcomes the distortions stemming from the undulations which ap- 
pear to be correlated across the different orders of structure functions (see Appendix B). 

Hints of the presence of log-periodic undulations decorating the power laws can be observed 
in Fig. 8 of Ref. [32] which studied the multifractal nature of turbulent energy dissipation, in 
Fig. 5.1 p.58 and Fig. 8.6 p.128 of Ref. [JT3J], Fig. 3.16 p.76 of Ref. [§], Fig. lb of Ref. g 
and Fig. 2b of Ref. ^ (to list a few). However, these barely perceptible undulations are far from 
proving the existence of a genuine discrete hierarchical cascade in three-dimensional flows, a la 
Richardson-Kolmogorov, as they are in general considered as "noise". Furthermore, if it exists, the 
scale ratio 7 has not be constrained. 

A preliminary attempt has dealt with another type of turbulence, namely freely decaying 2-D 
turbulence [26], in which experimental evidence and theoretical arguments suggest that the time- 
evolution of freely decaying 2-D turbulence is governed by a discrete time scale invariance rather 
than a continuous time scale invariance. The number of vortices, their radius and their separation 
display log-periodic oscillations as a function of time with an average log-frequency of ~ 4 — 5 
corresponding to a preferred scaling ratio of ~ 1.2 — 1.3. Since the statistical significance for each 
data set taken separately is not good when comparing with a large variety of different stochastic pro- 
cesses taken as null-hypothesis [53], Ref. [26] argued about a good statistical significance of this 
log-periodic oscillations on the basis of the coincidence of three highest peaks of the periodogram 
analysis which occur at log-frequencies far from the most probable frequencies occurring solely 
from noise. An alternative possibility is that the log-periodicity results from the long-range corre- 
lations in fractional Gaussian like-noise with Hurst exponent of the order of H 0.3. This case 
stresses the difficulty and complication of determination of the significance level of a log-periodic 
signal. 

The theory and practice of log-periodicity and its associated complex exponents has advanced 
significantly in the last few years [147], |3(J. It has become clear in the 1980s that complex expo- 
nents and log-periodicity appears as soon as a given physical problem is formulated on a discretely 
hierarchical geometry or network. More interestingly, only recently has it been realized that dis- 
crete scale invariance (DSI) and its associated complex exponents may appear "spontaneously" in 
Euclidean systems [41], i.e., without the need for a pre-existing geometrical hierarchy. Systems 
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where self-organized DSI has been reported include Laplacian growth models p3| , 17], rupture in 
heterogeneous systems [|], ^3[ ], earthquakes [[l4|, [if], |(], |T^], animals pT|], biased diffusion in 



percolating systems [pOfl as well as in financial time series preceding crashes H22L 24, 11, 49]. In 



addition, general field theoretical arguments | |4l| ] indicate that complex exponents are to be ex- 
pected generically for out-of-equilibrium and/or quenched disordered systems. This together with 
Novikov's argument suggest to revisit log-periodicity in turbulent signals. Demonstrating unam- 
biguously the presence of log-periodicity and thus of DSI in turbulent time-series would provide an 
important step towards a direct demonstration of the Richardson-Kolmogorov cascade or at least of 
its hierarchical imprint. 

The present paper presents the results of such an investigation on the log-periodic oscillations 
of the energy dissipation rate in 3-dimensional fully developed turbulence. In section 2, we explain 
where lie the difficulties of this enterprise and what are the classical traps of standard averaging 
methods. We stress in particular that standard ensemble averaging destroy the putative log-periodic 
undulations in the limit of large sample sizes. To address this problem, we explain our new "canon- 
ical" averaging methodology for extracting a possible log-periodic signal in turbulence time series 
and test it on synthetic examples. In section 3, we present the result of our analysis for different 
orders of the moments of the energy dissipation rate. In particular, we discuss the sensitivity of the 
results with respect to the filtering parameters of our technique. Appendix A presents a study of the 
dependence of the bandwidth selectivity of the Savitsky-Golay filter as a function of its parameters 
M (order of polynomial fit) and 2Nl + 1 (number of points) using the smoothing procedure. Ap- 
pendix B proposes a simple multifractal model accounting for the observation that the same set of 
log-frequencies, within numerical accuracy, are found independently of the moment orders varying 
from q = 1 to 19. 



2 "Canonical" averaging of undulations with noisy phases 

2.1 Characterization of the effect of averaging on the detection of log-periodicity 

If discrete scale invariance (DSI) exists in turbulence, this should be reflected in the existence of 
a hierarchy of characteristic scales £,£/j,£/'j 2 , ... where 7 is a preferred scaling ratio and £ is a 
macroscopic scale. 

Actually, this naive picture has to be revised to account for dynamical effects. If it exists, let 
us imagine that the energy cascade starts over some region at some instant from some scale £ a 
and then proceeds down the scales over a few generations of the hierarchy. It will probably be 
interrupted by various competing processes that break down the localness of the interactions due 
for instance to the influence of very thin shear layers [^] or slender vortex filaments [29], 80. Later, 



at another time in some other region, another energy cascade may start from some other scale £f, and 
then proceeds downwards over a different number of generations of the hierarchy. If this picture 
has some element of truth, there is an immediate consequence with respect to the detectability of 
log -periodicity associated with these transient cascades: 

1. such log-periodic undulations will be transient in a long time series; 

2. different burst of log-periodic oscillations will have different phases. 
Indeed, consider a pure log-periodic signal 

S(r) = cos(27r/ lnr — ip a ) (1) 
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of a structure function or moment as a function of scale r, where 



/ = l/ln 7 (2) 

is the log-frequency. The phase ip a can be redefined as ip a = 2nf In i a such that the undulation 
becomes cos(27r/ ln(r/^ a )). The values ro = £ a ,ri = £ a /lif2 = ^a/l 2 , • •• of the scale r corre- 
sponding to the maxima of the cosine function provide the discrete scales of the underlying cascade, 
which are all proportional to the "root" length scale £ a . If a second cascade is triggered from a 
slightly different scale the resulting oscillation cos(27r/ lnr — tpt,) with the same log-frequency 
will be out-of-phase compared to the first one with a phase shift ip a — tpb = 27r/ ln(£ a /£j,). This 
scenario suggests that, while the physics of the cascade may be universal with an universal scaling 
factor 7, the phases are sensitively dependent on the specific scale from which a given transient 
cascade nucleates. As the cascade may be dynamically triggered from a variety of scales £ a from 
the integral scale and deep in the inertial range, the corresponding phases are expected to be non- 
universal. Averaging a signal over several such transient cascades will thus wash out the information 
by "destructive interferences" of the log-periodic oscillations, giving them the appearance of noise. 

A similar situation has already been documented in other systems. Numerical simulations on 
Laplacian growth models [^, [lTj, 21] and renormalization group calculations [41] have taught 



us that the presence of noise or disorder modifies the phase in the log-periodic oscillations in a 
sample specific way leading to a "destructive interference" upon averaging. As sample-to-sample 
log-periodic undulations are destroyed by averaging over a long record, this is related to the self- 
averaging property |7],|3(], |52| ]. 

2.2 Existing methods for detecting log-periodicity 

Up to now, four different methods have attempted to improve on standard averaging techniques 
which, in our context, may "throw away the baby with the bath". 

1 . A first method consists in fitting data with a power law decorated by the log -periodic cor- 
rections as done to improve the determination of the time of occurrence of rupture [jl|, [23| ], 
of earthquakes [@, f«| |(J ||] and of financial crashes ^ [24], g|]. This parametric 



approach has the drawback of becoming unstable when too much noise is present. 

The second idea is to avoid performing any average altogether and, studying the main log- 
frequency of each realization, to construct their distribution over the whole set of realizations. 
This was used in an early attempt to demonstrate the presence of DSI structure in growing 



diffusion-limited-aggregation clusters [20]. 



The third non-parametric approach uses the fact that a periodogram obtained with the Lomb 
method [[57|] (see also below for more details) gives the power spectrum as a function of 
log-frequency independently of the phase of the undulations. As a consequence, if there is 
a genuine peak, even with a lot of noise, performing a Lomb periodogram for each sample 
and then averaging the Lomb periodograms will tend to make the genuine peak stand out of 
the noise if there are sufficiently many independent samples. This is the method used for the 



analysis of 2-D freely decaying turbulence Q26Q already mentioned. 



4. Ref. [21 ] have adapted to the problem of log -periodic oscillation the novel so-called "canoni- 



cal" averaging scheme proposed in [36] based on the determination of a realization-dependent 



effective critical point obtained from, e.g., a maximum susceptibility criterion. Conceptually, 



4 



the method re-sets the nucleus scale l a of the cascade in different realizations to approxi- 
mately the same value and then only afterwards perform the averaging. The non-parametric 
method has been tested successfully on diffusion limited aggregation clusters and on a model 



of rupture [|21||. For the problem of diffusion-limited aggregation, "canonical" averaging 
amounts to average at fixed numbers of particles per cluster. In the rupture problem, a max- 
imum rate of released energy is used to determine the position of the realization-specific 
critical control parameter which is then used to perform an average over many different real- 



izations at fixed reduced distance from the critical rupture. Ref. [ |48fl suggested to adapt this 
idea to turbulence by dividing a long turbulence time series into intervals of length equal to 
a turn-over time and rephasing the structure functions estimated over each such time interval 
by identifying a scale at which the dissipation rate is locally the largest. 

2.3 Our new improved "canonical" averaging procedure 

Here, we present a novel extension of the fourth method that turns out to perform extremely well. 
We now describe and illustrate it on simple synthetic data before turning to the full-Hedge analysis 
of a long turbulent time series. Since our goal is to detect log-periodicity in the structure function 
for which the natural variable is the logarithm of the scale, we can rephrase the problem into one of 
detecting a periodic undulation in a noisy time series. 

Consider for the sake of presentation, a signal similar to (|]). Posing t = lnr, we have the 
following "time" series 

y{t) = cos(27r/£ + ^(i)) , (3) 

where the phase ip(t) is now time varying to account for the variation of the scale of the nucleus 
of the cascade described above. Here, we do not show tests with additive noise as this standard 



situation was studied in [37, 18] and is relatively easily addressed with the Lomb periodogram 
spectral analysis (see also below). We thus focus our attention on the most difficult situation of 
random phases. 

To be specific, let us assume that ift(t) undergoes a random walk 

i/,(t) =ip(t-dt) + 2irs/J Amn(t) , (4) 

where ran(t) is a Gaussian random number of zero mean and of variance dt. With this parameteri- 
zation (|f), over N periods 1// of the deterministic component, the phase wanders by an amplitude 
of the order of 2-kA^/~N which is rapidly larger than 2ir if A is comparable to 1 as we discuss 
here. This random phase modulation can also be seen as a multiplicative noise acting on the sig- 
nal (cos(2tt ft), sm(2ir ft)) in the complex plane since cos(2ir ft + tp(t)) = cosV> cos(27r/t) — 
simp sin(27r/t). In the absence of the deterministic signal 2irft in the cosine in expression (||), 
the random process y(t) has a correlation function which can be evaluated exactly: (cos(t^(t + 
t)) cos(V'(t))) — (cos(ip(t + r))) (cos (?/>(£))) = e~ 2n T with a characteristic correlation time 
l/2ir 2 A 2 . 

Figure [j] shows one realization with 100 points of the time series ([3|,|]) with f = 1, A = 0.3 
generated numerically with dt = 0.1 in the time interval [0, 10]. This corresponds to a correlation 
time \/2ir 2 A 2 = 0.44. This value A = 0.3 thus corresponds to a very strong phase disorder, 
so strong that the time series appears very random. Figure [2| shows a standard spectral analysis 
(Fourier spectrum) of a long time series with 2000 points in the time interval [0, 200]. It is clear 
that the random walk of the phase ip(t) destroys completely the information of the existence of an 
underlying deterministic component at the frequency / = 1: the figure shows a completely noisy 
spectrum. 
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Figure |3] shows the result of the third method listed in section 2.2: we take the long time series 



with 2000 points in the time interval [0, 200] and divide it in 20 contiguous time series of 100 points 
each that are similar to that represented in Fig. |I[ We perform on each of these 20 time series a Lomb 
periodogram analysis and then average these 20 Lomb periodograms. The Lomb periodogram 
analysis performs local least-square fits of the data by sinusoids centered on each data point of the 



time series []42|]. The Lomb periodogram usually out-performs significantly other spectral methods 
for unevenly sampled and relatively short signals [ ]4"2| , |l6| , [37| , |53| ]. Figure || provides already a 
significant improvement over the standard Fourier analysis of Fig.g as a quite broad but rather clear 
peak centered approximately on / = 1 can be distinguished. However, its statistical significance 
is not convincing since the hypothesis that it could result just from random occurrences can not be 



rejected by standard statistical tests [ 37 , 53 ] 



Our new method consists in the following steps. 

1. Rather than analyzing the long time series of 2000 points, we divide it again in 20 contiguous 
time series of 100 points each: [1, 100]; [101,200]; ...[l+nxl00, (n+l)xl00]; ...; [1901,2000] 
These two numbers 20 and 100 are rather arbitrary and can be modified by and large. We 
need however a sufficient number of times series to average over as discussed below and a 
sufficient number of points in each time series to carry at least a few oscillations. 

2. On each sub-time series of 100 points, we apply the Savitsky-Golay filter [ [37| ] over a running 
window of size Nl+Nr+1 (with Nl points to the left and Nr points to the right of a running 
point) by fitting a polynomial of order M. We take Nl = Nr (symmetric filter) and vary Nl 
between 5 and 10 and M between 4 and 7, both by unit increments. This provides a total of 
24 filtered versions of each sub-time series of 100 points. The Savitsky-Golay filter has the 
advantage of providing smoothing without loss of resolution and of avoiding the distortion of 
the moments of the signal of all orders up to the order of the polynomial. We vary Nl and 
M over rather broad intervals to check for the sensitivity of this filtering process and to avoid 
the occurrence of spurious frequencies that may be created by the filter (see Appendix A). 

3. For each of the 24 filtered versions of each sub-time series of 100 points, we use its local 
polynomial representation to calculate analytically its local derivative by simply differentia- 
tion the polynomial associated with each point along the time series. The step constructs an 
observable which is sensitive to the phase. 

4. For each time series of the derivative, we identify the time i max of the local maximum closest 
to its middle point. This is analogous to identifying a local maximum of the control parameter 



of the "canonical" averaging method described above [Bq, 21] 



5. We then translate the origin of time of each derivative time series to this time t max - This is 
the "re-phasing" step. 

6. For a given Nl and M, we perform the average of the 20 time series derivatives. Varying Nl 
between 5 and 10 and M between 4 and 7, both by unit increments, provides a total of 24 
averaged time series derivatives that are shown in Fig. ||. This is the "canonical" averaging 
step. The figure is not unlike the local correlation structure g(r) obtained in the probing 
of the local order of liquids around a representative molecule. In particular, the oscillations 
decay in amplitude as a function of the distance to the central peak. Note that our choice of 
rephasing the local maximum closest to its middle point (and not other point) of each time 
series provides the best choice in order to observe a maximum number of oscillations. 
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7. For each of the 24 averaged time series derivatives, we perform a Lomb periodogram analysis 
in order to extract the most probable frequency of each signal. Figure |5] shows the correspond- 
ing 24 Lomb spectra. A significant fraction of the 24 Lomb spectra exhibit a very strong and 
statistically significant thin peak at the correct frequency / = 1 ± 0.05. Other spurious peaks 
are also present: they result from the interplay between Savitsky-Golay smoothing procedure 
and the phase noise. 

8. To confirm that the information on the existence of the deterministic signal at the frequency 
/ = 1 is present in each of the 24 Lomb spectra, Fig. ^| shows the average over these 24 
Lomb spectra in Fig. |5[ A very clear peak at / = 1.03 ± 0.1 is obtained with an amplitude 
twice as big as the background. 

We have presented here our new method on an extreme case of phase disorder. For smaller 
disorder, for instance already with A = 0.1, the signal obtained with our method is overwhelming 
while the direct Fourier spectral analysis or even the third method of direct Lomb averaging provide 
still unconclusive proof of the existence of a deterministic frequency. 

We now present the application of this technique for the detection of log-periodicity in the 
moments of the energy dissipation rate. 



3 Detection of log-periodicity in the moments of the energy dissipation 
rate 

3.1 Standard preliminary tests on the experimental data 

Very good quality high-Reynolds turbulence data have been collected at the S 1 ONERA wind tun- 
nel by the Grenoble group from LEGI [0]. We use the longitudinal velocity data obtained from 
this group to check for the possible presence of log-periodic undulation in moments of the energy 
dissipation rate. Figure [7] shows a typical sample of the time series of the streamwise component of 
the flow velocity, denoted as v hereafter in this work. 

The mean velocity of the flow is approximately (v) = 20ms~ 1 (compressive effects are thus 
negligible). The root-mean-square velocity fluctuations is v Tms = l^ms^ 1 , leading to a turbulence 
intensity equal to / = v rms / (v) = 0.0826. This is sufficiently small to allow for the use of Taylor's 
frozen flow hypothesis. The integral scale is approximately 4m but is difficult to estimate precisely 
as the turbulent flow is neither isotropic nor homogeneous at these large scales. 



The Kolmogorov microscale r\ is given by [[32|] rj 



u 2 {v) 2 



1/4 



0.195mm, where v 



15 {(dv/ dt) 2 ) _ 

1.5 x 10 _5 m 2 s _1 is the kinematic viscosity of air. dv/dt is evaluated by its discrete approximation 
with a time step increment dt = 3.5466 x 10~ 5 s corresponding to the spatial resolution bi = 0.72mm 
divided by (v). 



The Taylor scale is given by [32] A = 77^7^^172 = 16.6mm. The Taylor scale is thus about 



85 times the Kolmogorov scale. The Taylor-scale Reynolds number is Re\ = VzmA = 2000. This 
number is actually not constant along the whole data set and fluctuates by about 20%. 

We have analyzed a total number of about M = 1.4 x 10 7 data points provided in 200 segments 
of 2 16 = 65536 data points. In order to implement our detection method presented in section 2.3, 
we have re-partitioned the total data set of N points into N "records" of length L. In the sequel, we 
present results for N = 100 records of length L = 2 17 130, 000 data points (N x L = AO and for 
N = 20 records of L = 5 x 2 17 ps 655, 000 data points. Having N = 100 (or N = 20 respectively) 
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such records will allow us to test for the reproductivity and stationarity of the results. Each of the 
N records of length L is then itself subdivided into "samples" of 2 11 = 2, 048 data points. We thus 
have 64 (respectively 320) samples per record for L = 2 17 (respectively for L = 5 x 2 17 ). Taking 
an integral scale of £q = 4m with one point per 0.72mm, a complete turnover covers approximately 
21q/8i 10 5 data points. Each sample thus covers a finite fraction (« 1/5) of a turn-over time. 

We have checked that the turbulent velocity time series recovers the standard scaling laws pre- 
viously reported in the literature with similar quality. In particular, we have verified the validity of 
the power-law scaling E(k) ~ with an exponent (3 very close to | over a range more than two 
decades, similar to Fig. 5.4 of [13] provided by Y. Gagne and M. Marchand on a similar data set 
from the same experimental group. Similarly, we have checked carefully the determination of the 
inertial range by combining the scaling ranges of several velocity structure functions Jl"4] ] (see also 
Fig. 8.6 of JT3|]). Conservatively, we are led to a well-defined inertial range 60 < £/r) < 2000. 



3.2 Moments of the energy dissipation rate 

Using Taylor's hypothesis which replaces a spatial variation of the fluid velocity by a temporal 
variation measured at a fixed location, the rate of kinetic energy dissipation at position i is 

ei ~ [(v i+ i - Vi) /5i} 2 , (5) 

where 8$. is the resolution (translated in spatial scale) of the measurements. The total dissipation 
rate Ei in a spatial domain Q of size £ is 

E t = 8 t e * > (9 

where fij is the ith sub-piece of linear dimension 5g in Q, such that n Qj ■ = $ for i ^ j and 
Uj Oj = f2. We use normalized energy dissipation Ee/ 'El, where L is the size of the system, and 
normalized scale £/rj. 

We study the q-th moment of the normalized energy dissipation rate 

M q (£) = J2(E e /E L y , (7) 

as a function of the scale £. The summation in (^) is performed over all domains of size I within the 
system of size L. The signature of multifractal scaling is that the scaling law 

M q (£) ~ (£/fl) riq) , (8) 
is expected to hold in the inertial range with moment exponents r(q) that are non-linear functions of 



q. Using standard arguments [B5L 15, 13], the multifractal spectrum /(a) is the Legendre transform 



of r(q). The validity of these power laws (||) and multifractality have been verified with a good 
accuracy in experiments [32, 13]. We have reproduced these results reliably. As already mentioned, 
one can discern sometimes some small noisy oscillations, which are suggestive of log-periodicity. 



3.3 Canonical averaging of the moment exponents 

In order to test for the possible existence of DSI, we apply the methodology explained in section 
[2~3| to each of the ./V records of length L = 2 17 . We estimate the moments M q {£) defined by (||) 
in the range £jr\ G [216, 1840] which is well within the inertial range for each sample. The lower 
boundary 216 is 3.6 times larger than the lower bound 60 of the inertial range previously discussed, 
in order to avoid the effect of very small scales that make the signal more noisy. 
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1. Given a record of length L = 2 17 (resp. L = 5 x 2 17 ), we cut it into 64 (resp. 320) samples 
of size 2 11 points. 

2. We then calculate the local moment exponent r (q, ln(£/r])) as the logarithmic derivative of 
M q (£) with respect to £jr\\ 

3. For each sample, we then identify the sample-dependent reduced scale £ c /rj closest to the 
central point of the scale-interval of the analysis at which r (q, ]n(£/rf)) is maximum. 

4. We rephase all functions r (q,ha.(£/r])) of the variable ln(£/r]) of each sample so that their 
£ c /r] coincide. This allows us to define the reduced variable A = \n.{£/rj) — ln(£ c /rj) such 
that all functions r (g, A) are rephased at A = 0. 

5. We then average these rephased functions r (q, A). This defines the function (local average 
exponent) 

D(A) = (r(q,A)) , (10) 

where the average is performed over the different samples in one record at fixed A = 
\a.{£/rj) — \n(£ c /rj) (where \n(£ c /rj) vary from sample to sample and ]n(£/rf) is adjusted 
accordingly). 

6. We then perform several comparative analysis to check the quality of the results discussed in 
the following. 

3.4 Results 

First, let us stress that our results presented below have been found to yield the same set of log- 
frequencies, within numerical accuracy, independently of the moment orders varying from q = 1 
to 19. Appendix B proposes a simple multifractal model explaining this observation. The standard 
deviation of the log-frequencies discussed below determined from the N = 100 records of length 
L = 2 17 shows that it is larger for q = 1 and reaches a plateau at q = 3. We have thus chosen to 
present the following results for q = 3. 

Figure || shows the canonically averaged local log-derivative D(A) defined in (|10|) as a function 
of A (left panels) and their corresponding Lomb periodograms (right panels) for three choices of 
the Savitsky-Golay filter parameters Nl and M. This analysis is performed over the 64 samples of 
a single record of length L = 2 17 points. 

Rather strong periodic oscillations can be observed on the local average exponent D(A), which 



results in the high Lomb peaks in the right panels which have good statistical significance [ p7 , 53]. 
As expected, the central oscillation at A = is the largest and the undulations on both side decay 
as a result of the progressively increasing destructive interference due to noise. One can also note a 
slight left-right asymmetry that we attribute to the stronger noise affecting the small scales (A < 0) 
compared to the large scales (A > 0). 

While the undulations and their Lomb spectra seem convincing, there is a problem that seems 
at first sight to destroy the evidence of a genuine log -periodicity: the periods of the oscillations 
are strongly dependent upon the parameters Nl and M of the Savitsky-Golay filter (see Appendix 
A): it would thus seem that we are quantifying spurious oscillations created by the massaging of 
the data. A similar effect has been found earlier associated with the construction of cumulative 
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quantities [|T8|]. In order to understand the possible origin of the strong variations in significant log- 
frequencies as a function of Nl and M and following Ref. [fl5|], we have collected the two most 
significant log-frequencies (two highest Lomb peaks) of all the 2400 Lomb periodograms available 
from our analysis of 100 independent records by the 24 possible filters (spanning Nl = 5 — 10 
and M = 4 — 7) used for each record. Fig. ||(a) (respectively (b)) shows the histogram of the 
log-frequencies associated with the largest (respectively second largest) Lomb peaks. The vertical 
dashed lines correspond to log-frequencies equal respectively to fa = 1.44, fa = 2/i = 2.89, fa = 
3/i = 4.33, /4 = 4/i = 5.77 and fa = 5/i = 7.21. These values correspond respectively 
to increasing harmonics of a fundamental frequency fa = 1/ In 7 associated with the scale ratio 
7 = 2. Remarkably, similar histograms (but with smaller amplitudes for the peaks at fa and f%) are 
also obtained by collecting the frequencies corresponding to the third and fourth largest Lomb peaks 
among the 2400 available Lomb periodograms. The existence of these 5 harmonic frequencies is 
thus a very robust feature. Another evidence of the quality of the description in terms of 5 harmonics 
of a fundamental frequency, we represent in Fig. 10 the value of the frequency measured at the nth 
maximum of the histogram shown in Fig. ^(a) as a function of the order n of the maximum. The 
straight line shows the excellent linear fit to the equation f n = nfa, where the sole adjustable 
parameter fa is found equal to / = 1.62 ± 0.1. This value corresponds to a prefered scaling ratio 
of the cascade equal to 7 = 1.85 ± 0.1. 

We are now in position for explaining the strong variations in the position of the Lomb peaks 
observed in Fig. ||. Since each filter has a different passing-band (see Appendix A), for a periodic 
signal which has harmonics of a fundamental frequency with strong amplitudes, these different har- 
monics will be differently expressed by the different filters. This is what we observe in Fig. ||. The 
two top panels correspond to the filter with (Nl = 10, M = 4) which selects preferentially the 
fundamental frequency fa. The two middle panels correspond to the filter with (Nl = 9, M = 5) 
which selects preferentially the second harmonic fa = 2/i. Finally, the two bottom panels corre- 
spond to the filter with (Nl = 8, M = 7) which selects preferentially the third harmonics fa = 3 fa. 
Furthermore, the fact that the second or third harmonics have a larger amplitude than the fundamen- 
tal frequency is not a surprise and is found for instance in multifractal measures constructed on DSI 
fractal geometries [64|]. 

The limited bandwidth of each filter and their selection of different harmonics is further quan- 
tified by Table [l] which explores the range of values Nl from 5 to 10 and M from 4 to 7, corre- 
sponding to a total of 24 filters. Each filter is applied to each of the N = 100 records of length 
L = 2 17 to obtain the canonically averaged D(A) from which we obtain the corresponding 100 
log-frequencies as explained in section 3.3. The average (/) of these 100 log-frequencies and their 
standard deviation Of are given as a function of M and Nl in Table [I]. (/) increases with increasing 
M and/or decreasing Nl, an expected property of the Savitsky-Golay filter (see Appendix A). This 
is not however where lies the potentially useful information. To retrieve it, we identify ten pairs 
(M, Nl), specifically (4,6), (4,7), (4, 10), (5,8), (5,9), (5,10), (6,8), (6,9), (6,10), and (7,8), 
whose cif < 0.3, i.e., such that there is good confidence in the determination of a log-frequency. 
Interestingly, these ten pairs give three groups of log-frequencies: (1) eight pairs marked by astroids 
(*) have an averaged log-frequencies in the range 2.71 — 2.99; (2) one pair marked by the star (*) 
has (/) = 4.37; (3) one pair marked by the circle (o) has (/) = 1.79. These three frequency 
bands are in good agreement with the values for the harmonics fa = 1.44, fa = 2 fa = 2.89 and 
/3 = 3/i = 4.33 already discussed. 

Note also that there are pairs, say (M = 5,Nl = 7) and (M = 6,Nl = 7), leading to 
very large standard deviations and values for (/) that are far from any of the frequencies fa, 
2/i,3/i,4/i,.... This apparent negative result actually hides a remarkable confirmation of our 
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N L = 5 


N L = 6 


N L = 7 


N L = 8 


N L = 9 


Nl = 10 


M = 4 


3.63 ±0.91 


2.89 ±0.20* 


2.71 ±0.27* 


2.51 ±0.49 


2.18 ±0.81 


1.79 ±0.20° 


M = 5 


5.22 ±0.63 


4.32 ±0.36 


3.76 ± 0.70 


2.99 ±0.25* 


2.85 pmO.U* 


2.76 ±0.16* 


M = 6 


5.22 ±0.64 


4.31 ±0.36 


3.78 ±0.69 


3.01 ±0.28* 


2.85 pmO.U* 


2.77 ±0.11* 


M = 7 


6.15 ±0.66 


5.72 ±0.63 


4.82 ±0.55 


4.37 ±0.27* 


3.. .90 ±0.60 


3.13 ±0.43 



Table 1: Averaged logarithmic frequencies (/) as a function of M from 4 to 7 and Nl from 5 to 
10 of the Savitsky-Golay filter. We use 6 x 4 = 24 versions of the Savitsky-Golay filter to test for 
the sensitivity of the order M of the fitting polynomial and the number 2Nl + 1 of points. This 
24 filters allow us to sweep broadly a large range of possible log-frequencies and thus to access the 
fundamental log-frequency f\ and several of its harmonics. 



*//</> 


N L = 5 


N L = 6 


N L = 7 


Nl = 8 


Nl = 9 


N L = 10 


M = 4 


0.25 


0.07* 


0.10* 


0.20 


0.37 


0.11 


M = 5 


0.12 


0.08* 


0.19 


0.08* 


0.05* 


0.06* 


M = 6 


0.12 


0.08* 


0.18 


0.09* 


0.05* 


0.04* 


M = 7 


0.11 


0.11 


0.11 


0.06* 


0.15 


0.14 



Table 2: Relative deviations a y / (/) of the averaged log-frequencies already shown in Table [I] as a 
function of the parameters M and Nl of the Savitsky-Golay filter. 



previous presentation. Indeed, Fig. 11 [a) shows the detected log-frequencies / of the 100 dif- 



ferent records calculated for different orders q of the moments of the energy dissipation rate for 
(M = 5, Nl = 10). Note first that the determined frequency for q = 1 fluctuates more than for 
higher moment orders q > 1, which justifies our investigation of different moment orders. Sec- 
ondly, the arrow shows a large fluctuation occurring for record 74: while most of the records are 
tuned to the frequency f'2 = 2/i, record 74 selects f± as its best log-frequency. This provides a 
mechanism in which the competition between the different first frequencies can lead to an average 
in between as found in Table [I] for (M = 5, N L = 7) and (M = 6,N L = 7). Figure |ll](b) shows 
all the Lomb periodograms obtained for record 74: we directly visualize in this case the strong com- 
petition between the two first log-frequencies f\ and /2, strengthening the case for the existence of 
a discrete set of log-frequencies. 

Table H| offers another perspective to this analysis by showing the relative deviations <Jf/(f) 
of the average log-frequencies shown in Table ^[ To illustrate the robustness of the recovery of 
the main harmonics of the fundamental frequency fx, let us select only those log-frequencies de- 
termined with a relative precision better than 0.10. We find now two classes of pairs (M, Nl) of 
filters. The first class contains (4,6), (4,7), (5,8), (5,9), (5, 10), (6,8), (6,9) and (6, 10) and is 
marked by the astroid (*) and corresponds to (/) ~ 2.71 — 2.99. The second class contains (2) 
(5, 6), (6, 6) and (7, 8) and is marked by the stars (*) and corresponds to (/) ~ 4.31 — 4.37. 

Increasing the size L of the records improves the quality of the signal. This is first quantified by 
Figs. 12-15 in which the dashed lines corresponding to I = 5 x 2 17 give larger and thinner peaks 
that the continuous lines obtained for L = 2 17 (see below). This confirms the role of our canonical 
averaging procedure which improves in quality as the number of samples used in the averaging 
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</>(<7/</» 


^L = 5 


N L = 6 


N L = 7 


N L = 8 


N L = 9 


Nl = 10 


M = 4 


3.68(0.20) 


2.90(0.03) 


2.82(0.04) 


2.67(0.08) 


1.92(0.18) 


1.78(0.05 ) 


M = 5 


5.38(0.11) 


4.38(0.02) 


4.07(0.12) 


2.96(0.03) 


2.84(0.02) 


2.82(0.02 ) 


M = 6 


5.53(0.08) 


4.38(0.06) 


4.06(0.12) 


2.96(0.04) 


2.84(0.03) 


2.80(0.02 ) 


M = 7 


6.34(0.11) 


5.74(0.02) 


4.72(0.10) 


4.40(0.02) 


4.33(0.02) 


2.90(0.04 ) 



Table 3: Log-frequencies (/) and relative deviations <Jf/(f) as a function of M and Nl for the 
records of length L = 5 x 2 17 . 



increases. Table || gives the log-frequencies (/) and their corresponding relative deviations 0"//(/) 
as a function of M and Nl of the Savitsky-Golay filter for the N = 20 records of length L = 5 x 2 17 . 
Comparing Table | (L = 5 X 2 17 ) with Tables [T] and | (L = 2 17 ), we find that a larger L leads to a 
clearer determination of the harmonics log-frequencies. For instance, keeping all frequencies such 
that the relative error (Jf/{f) is less than 0.1, we find that all the determined log-frequencies now 
very close the exact values of the harmonics of f c = 1.44. This confirms an important expectation 
of our canonical averaging procedure: using more samples will definitely improve its efficiency for 
extracting some hidden fundamental log-frequency and/or it harmonics. 

Figures summarize our analysis by presenting the Lomb periodogram obtained by averag- 



ing over the N = 100 (respectively N = 20) records of length L = 2 17 (respectively L = 5 x 2 17 ), 
for each of the 24 Savitsky-Golay filters described above. All the five harmonics of the fundamental 
frequency fi discussed above are seen to correspond to the main peaks of the Lomb periodograms 
with a good approximation. It is expected that their amplitude should vary with the filter, in re- 
sponse to the frequency selection of each filter. It is quite remarkable however to obtain such a 
strong coherence of the peaks over all the filters. This coherence is further demonstrated by taking 
the average over all 24 x N Lomb periodograms (24 filters for each of the N records) which is 
shown in Fig. [T^. Notice the good agreement between the values of the log-frequencies corre- 
sponding to the peaks and the prediction from DSI with a root scaling ratio 7 = 2 (shown as the 
dashed vertical lines). We attribute the only significant deviation observed for the first harmonics 
fi to a combination of two factors: (1) the amplitude of the Lomb peak is small and its maximum 
is thus mechanically biased by the existence of the second large peak at $2', (2) as we have seen, 
different filters have different bandwidth and the averaging of a relatively low signal will lead to 
unavoidable distortions. 

3.5 Statistical significance 

The statistical significance of a Lomb peak requires in principle the specification of the properties 
(distribution function and dependence) of the noise decorating the signal. The simplest and most 
often used noise is the Gaussian white noise for which the statistical significance level of the highest 



Lomb peak can be determined analytically Q42| , |16| , |37[ ]. It gives an approximately exponentially 
decaying false-alarm probability: 

Pr{> z) = I - (1 - e' z ) M ' ~M'e~ z , (11) 

where M' is proportional to the number of data points. 

In practice, noise is almost never Gaussian white noise with often much fatter tails of the dis- 
tribution and possible short or long-range correlations. Extensive numerical simulations have in- 
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vestigated the impact of heavy-tailness, of correlations and of the interplay between them on the 
significance level of a Lomb periodogram peak [53]. It was found that the worst case leading to 
high false-alarm probabilities correspond to a persistent fractional Brownian noise with large Hurst 
exponent. Consider for instance the highest Lomb peak with height 16.5 presented in Fig. 14 [f) as 
the dashed line, obtained with 52 data points. According to the simulations of Ref. [53], this can 
be translated into a false-alarm probability of 0.1% and a statistical confidence of 99.9%, ever for 
the worst case of a strongly persistent fractional Gaussian noise with Hurst exponent H = 0.9. For 
a Lomb peak of height 10 as shown in Fig. [f4](b) as the dashed lines, the confidence level is about 
99.95% for a fractional Brownian noise with H = 0.7 and ~ 0.99 with H = 0.8. 

We have performed additional simulations to test whether the five log-frequency harmonics ob- 
tained from the Lomb peaks in Figs. 12-15 of the turbulent data could result from the strongest case 
of spurious log -periodicity [53, 18], namely when the noise exhibits a strong persistence with long- 
range correlations. We would like to compare the probability for obtaining a large Lomb peak in a 
purely synthetic fractional Brownian noise with a large Hurst exponent, using exactly the canonical 
averaging procedure used for analyzing the turbulence data. It is known [|3|, |l8|] that, the larger H 
is, the stronger are the spurious log-periodic undulations. Specifically, we use fractional Gaussian 
noises of Hurst exponent H = 0.99 to replace the original fluctuations decorating the leading power 
law of the moments for each sample and follow exactly the same analysis procedure as for the tur- 
bulence data. We generated 100 data sets of fractional Gaussian noises with the Hurst exponent 
H = 0.99, each set having 61 data points. The Lomb periodograms averaged over 100 samples, 
corresponding to different M and Nl of the Savitsky-Golay filter, are shown in Fig. 17. Compared 
with Figs. 12-15, the main feature of the Lomb periodograms in Fig. 17 is that there are no har- 
monics visible and the extracted frequencies appear completely random. This is in strong contrast 
with the analysis for the turbulence data which exalts the five log-periodic frequency harmonics. 
Moreover, the histogram of the total 2400 log-frequencies corresponding to the highest peaks of the 
24 x 100 Lomb periodograms, shown in Fig. |8|, does not indicate that any frequency is playing a 
special role. The same qualitative results are found when varying the Hurst exponent H between 
0.5 to 0.99. 

In summary, while the Savitsky-Golay filter leads to strong spurious log-periodic components 
in the fractional Gaussian noises with high Hurst exponent, their large dispersion from one filter to 
the other rules out that such a noise can be at the origin of the log-periodicity uncovered from our 
analysis. In addition, since the fractional Gaussian noises with high Hurst exponents are by far the 
processes leading to the largest spurious log-periodicity, this test also excludes all reasonable types 
of noise investigated previously [53, 18], such as for GARCH(1,1) noise with Student distributions 
with different numbers of degrees of freedom, for power law noises, for Levy stable noise as well 
as of course for independent Gaussian noise. 



4 Conclusions 

We have presented the results of a search for the existence of log-periodic oscillations decorating 
the moments of the energy dissipation rate in experimental data at high Reynolds number 2500) 
of three-dimensional fully developed turbulence. We have proposed a simple dynamical represen- 
tation of the Richardson-Kolmogorov cartoon of a cascade that explains why standard averaging 
techniques erase by their construction the possible existence of log-periodic corrections to scaling 
associated with a discrete hierarchy. We have introduced a novel "canonical" averaging method- 
ology for extracting a possible log-periodic signal in turbulence time series and have tested it on 
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synthetic examples constructed to capture the interplay between a weak log-periodic component 
and rather strong multiplicative and phase noises. The many tests presented confirm the remarkable 
observation of statistically significant log-periodic correction to scaling, with a prefered scaling 
ratio for length scales compatible with the value 7 = 2. A strong confirmation of this result is 
provided by the identification of up to 5 harmonics of the fundamental log-periodic undulations, 
associated with up to 5 levels of the underlying hierarchical dynamical structure. 
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provided by Y. Gagne. We are grateful to J. Delour and J.-F. Muzy for help in pre-processing these 
data. This work was partially supported by NSF-DMR99-71475 and the James S. Mc Donnell 
Foundation 21st century scientist award/studying complex system. 
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(f)/fc 


A^L = 5 


N L = 6 


N L = 7 


N L = 8 


N L = 9 


Nl = 10 


M = 4 


1.87 


1.82 


0.81 


1.02 


1.00 


1.04 


M = 5 


3.04 


1.87 


2.04 


1.82 


1.07 


0.98 


M = 6 


3.04 


2.04 


2.00 


2.00 


1.87 


1.42 


M = 7 


3.11 


3.02 


1.99 


1.95 


2.04 


1.99 



Table 4: Ratios of the frequencies (divided by / c ) extracted from the largest peak of the average 
of 10 Lomb periodograms of a synthetic signal with 4 harmonics and a large centered Gaussian 
white noise with standard deviation a e = 3 as a function of the two parameters M and Nl of the 
Savitsky-Golay filter. 

Appendix A: Study of the bandwidth selectivity of the Savitsky-Golay 
filter 

By construction, the Savitsky-Golay filter has a finite bandwidth, i.e., it not only filters out high 
frequencies but also enhances preferentially an intermediate frequency band whose determination 
depends on its two parameters, the order M of the fitting polynomial and the number 2Nl + 1 over 
which this fit is performed. Let us make this remark quantitative. 

Consider a periodic function with four harmonics of the fundamental frequency f c 

y(t) = cos(27r f c t) + cos(47r/ c t) + cos(67r/ c t) + cos(8ir f c t) + e, (12) 

where f c = l/ln2, t is evenly spaced in [—1, 1] with 61 points and the noise e is normally dis- 
tributed with zero mean and standard deviation a e = 3. This choice is reasonably close to the 
situation found the turbulence data, in that several harmonics and noise compete with similar am- 
plitudes. 

We generate 10 times series of 61 points with t regularly spaced in 6 [—1,1]. For each of these 
times series, we choose one pair (M, Nl) and apply the Savitsky-Golay filter. We then apply the 
Lomb periodogram technique on the resulting smoothed time series to obtain a Lomb spectrum. We 
average the 10 Lomb spectra and identify the frequency of the largest peak which is reported in the 
Table f|. Actually, we show the ratios of the extracted frequencies to f c . Remarkably, 17 frequencies 
among the 24 (i.e., 71%) are found to within less than 10% from one of the harmonics f c , 2/ c , 3/ c . 
The 7 remaining frequencies are within less than 20% from one of the harmonics / c , 2/ c , 3f c . The 
second observation is that, indeed, a given filter (M,Nl) amplifies selectively on average one 
harmonics at the expanse of the other. Interestingly, when increasing the noise amplitude a e , the 
fourth harmonic 4/ c is also sampled by a few filters. Obviously, when the amplitude of the noise is 
very large, spurious frequencies are selected but come with a low statistical significance. With large 
noise and some combinations of the amplitudes of the harmonics in ([Of), the artifactual frequency 
f m = 1.5/ In 2 appears; its origin has been unravelled in Ref. |19|] as corresponding to the most 
probable noise structure. For very large noise, we can also identify a kind of effective non-linear 
coupling between the frequencies leading to combinations of the harmonics 

We have performed many tests with pure noise, with a single cosine plus noise and with several 
harmonics plus noise. In the first case, the frequencies obtained as the maxima of the Lomb spectra 
do not exhibit any prefered value and fluctuate broadly from realization to realization. With a single 
oscillatory signal, we see a clear peak in the Lomb spectrum, broadened by the disorder. In the 
case of several harmonics of a fundamental frequency, we find that one of the harmonics in general 
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dominates each Savitsky-Golay filter, allowing us to detect all the relevant harmonics by using many 
different filters. 
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Appendix B: A simple multifractal model with discrete scale invariance 
(DSI) independent of the moment order q 

Our analysis of the log-periodic undulations decorating the power law of the moments of the dis- 
sipation rate suggests that the log-frequency of these undulations is independent of the order of 
the moments used in the analysis. Here, we propose a simple multifractal model with DSI which 
exhibits this property. 

Consider a simple cascade process of the type introduced by the Russian school and studied by 



Mandelbrot (see [ |12| ] and references therein), with however one simple modification emphasizing 
DSI. Each large eddy T of size r is assumed to transfer the fraction rrii of its energy to to n small 
eddies T\ of size Vi = r/%, i = 1, 2, • ■ ■ , n, where the pair (m^, 7j) each are one of n different fixed 
values (mi, 71), (m2, 72), (rn n ,j n ) respectively. Each small eddy of size r{ in turn transfers a 
fraction rrij to each of its n daughters of size r.j/77, where rrij and 7j are drawn from the same fixed 
set (mi, 71), (m-2, 72), (m n , j n ). This deterministic multiplicative cascade process continues 
down to the dissipation scale and creates a self-similar multifractal measure with multipliers and 
characteristic scaling ratios rrii and 7$, respectively. We do not restrict the model to obey the space- 



filling condition Y2=l l hi = 1 & 



Using the self-similarity of the cascade process, we can write 

n 

M q {F;£)=J2M q (fi;£) , (13) 



i=l 



where M q (tF; I) is the gth-order moment of the energy dissipation rate defined on eddy T . Accord- 
ing to the self-similarity of the multinomial measure, we have 

M q (F ) ;£) = m q i M q (? i ; li £). (14) 



Combining Eqs. ( 13) and (N4) and removing the reference to a specific eddy T gives the following 



discrete scale invariant equation [54] 



M q {l) = Y J M q { li £)m\ . (15) 



i=l 



For statistically self-similar measures, Eq. ( |15| ) is modified by inserting in each term at the r.h.s. of 
equation (|l5j) the specific probabilities of the realizations among the n scenarios at each step of the 
cascade, 



Equation (15) generalizes the equation M q (£) = m q M q (^£) embodying the discrete scale in 



variance of simple fractals constructed using the iteration of a deterministic rule such as the triadic 



Cantor set or the Sierpinsky gasket. This equation is discussed in Ref. [3C] in relation with so-called 
"fractal strings". Looking for a scaling solution M q {£) oc £ T ( g \ the exponents r(q) are the solutions 
of 

n 

E^ = 1 (16) 

1=1 

and form an infinite countable set of complex numbers. For the "lattice" case where there is a 
common root 7 and there are integers k{ such that 

li = 7 fcl , (17) 



the solution of ( |15| ) can be explicited analytically. 
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Substituting Eq. (17) into Eq. (|I5|), we have 

n 

M q (£) = J2M q (^£)mf, (18) 
»=i 

whose solution is 

M q (£) = C(£)£ T ^ , (19) 

where C{£) is a universal function and the case of C{£) = const, recovers the familiar power law 
given in Eq. (Q). It is easy to verify that if C{£) is a periodic function of the variable ln(£) with the 
log-period ln(7), i.e., 

C(£) = C{~i k H) , (20) 

then Eq. ([H]) is the solution of Eq. (|15|). Without loss of generality, we can rewrite the general 
solution as 

M q {£)=iP q M£)/ln{ 7 )]£ r M . (21) 

This implies that M q {£) is log-periodic. We emphasize that ijj q is (/-dependent, while the log -period 
ln(7) of M q {£) is independent of q. Another important implication is that, given a unique prefered 
scaling ratio 7 = exp(l//), it is possible to have a set of characteristic scaling ratios 7^ satisfying 



the lattice condition expressed by Eq. (17) which describe an infinite recursion of the multifractal 
measure. 
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Figure 1: One realization with 100 points of the time series with f = 1, A = 0.3 generated 
numerically with dt = 0.1 in the time interval [0, 10]. 
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Figure 2: Standard spectral analysis (Fourier spectrum) of a long time series ([||]) with / 
0.3 generated numerically with dt = 0.1 in the time interval [0, 200] with 2000 points. 
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Figure 3: Result of the third method listed in section 23. applied to the time series of 2000 points 
generating to obtain the spectrum shown in Fig. |2| We divide the long time series with 2000 points 
in 20 contiguous time series of 100 points each that are similar to that represented in Fig. [|. We 
perform on each of these 20 time series a Lomb periodogram analysis and then average these 20 
Lomb periodograms. 
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Figure 4: For a given Nl and M, we perform the average of the 20 time series derivatives described 
in the text. Varying Nl between 5 and 10 and M between 4 and 7, both by unit increments, provides 
a total of 24 averaged time series derivatives that are shown here. 
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Figure 5: For each of the 24 averaged time series derivatives represented in Fig. ^ we show the 
corresponding Lomb periodogram spectrum. 
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Figure 6: Average over the 24 Lomb spectra shown in Fig. |[ 
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Figure 7: (a) A typical streamwise velocity time series recorded by hot-wire anemometry in the SI 
ONERA wind tunnel, (b) Corresponding time series of the velocity increments. 
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Figure 8: Local log-derivative D(A) defined in (JlO|), canonically averaged over the 64 samples of 
a single record of length L = 2 17 points, as a function of A (left panels) and their corresponding 
Lomb periodograms (right panels) for three choices of the Savitsky-Golay filter parameters Nl and 
M. (a-b) M = 4 and N L = 10, (c-d) M = 5 and N L = 9, (e-f) M = 7 and N L = 8. The dotted 
line in the left panels are pure cosine functions D(A) = (D) + \f2ao cos(27r/A). 
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Figure 9: Histograms of all log-frequencies corresponding to (a) the most significant and (b) the 
second most significant Lomb peaks of all 2400 Lomb periodograms available from our analysis 
of 100 independent records by the 24 possible filters (spanning Nl = 5 — 10 and M = 4 — 7) 
used for each record. The vertical dashed lines correspond to log-frequencies equal respectively 
to h = 1.44, f 2 = 2/i = 2.89, f z = 3h = 4.33, U = 4/i = 5.77 and f 5 = 5 A = 7.21, 
which correspond to increasing integer powers of a fundamental scale ratio exactly equal to 7 = 2, 
through the relationship f\ = 1/ In 7. 
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Figure 10: Value of the frequency measured at the nth maximum of the histogram shown in Fig. 
| as a function of the order n of the maximum. The straight line is a linear fit with the equation 
fn = nfi, where the adjustable parameter /i is found equal to fx = 1.62 ± 0.1. 
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Figure 11: (a) Leading log-frequencies / of different records (numbered on the abscissa) calculated 
for different moment order q marked with different symbols (q = 1 corresponds to the open circles). 
The arrow indicates a large fluctuations of the best log-frequency occurring for the 74th record, (b) 
All the Lomb periodograms obtained for record 74. 
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Figure 12: Lomb periodogram obtained by averaging over the iV = 100 (solid curve), respectively 
the N = 20 (dashed line), records of length L = 2 17 , respectively L = 5 x 2 17 . The order of the 
fitting polynomial of the Savitsky-Golay filter is fixed at M = 4, while iV^ goes from 5 to 10 in 
panels (a) to (f). 
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Figure 13: Same as Fig. 12 with the order of the fitting polynomial of the Savitsky-Golay filter 
fixed at M = 5. 
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Figure 14: Same as Fig. 12 with the order of the fitting polynomial of the Savitsky-Golay filter 
fixed at M = 6. 
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Figure 15: Same as Fig. 12 with the order of the fitting polynomial of the Savitsky-Golay filter 
fixed at M = 7. 
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Figure 16: Average of all 2400 (continuous line) (resp. 480 (dashed line)) Lomb periodograms 
shown in Figs. 12-T3] over all 24 filters and all records. The vertical dashed lines correspond to log- 
frequencies equal respectively to f\ = 1.44, f'2 = 2/i = 2.89, fa = 3/i = 4.33 and f± = 4/1 = 
5.77. These values correspond respectively to increasing harmonics of a fundamental frequency 
fi = 1/ In 7 associated with the scale ratio 7 = 2. 
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Figure 17: Lomb periodograms averaged over 100 samples of fractional Gaussian noises with the 
Hurst exponent H = 0.99, each set having 61 data points. One average Lomb periodogram is 
shown for each of the 24 values of the pair of parameters of the Savitsky-Golay filter: the order M 
of the fitting polynomial of the Savitsky-Golay filter ranges from 4 to 7, while Nl varies from 5 to 



10. These average Lomb periodograms should be compared with Figs. 12-15 
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Figure 18: Histograms of all log-frequencies extracted from the largest peaks of the 2400 Lomb 
periodograms of synthesized fractional Gaussian noises with H = 0.99. 
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